function int_cost_tax = fn_int_cost_app_tax(R1,delta,D0,chi,~,prob,phi,dist_s,s_min,~,dist_d,d_min,d_max,s_hat,s_star,~,kappa)
% This function calculates individual welfare

f_int_tax = @(s) fn_dcdd_tax(delta,R1,D0,chi,s,phi,dist_d,d_min,d_max,kappa).*dist_s.pdf(s);

A_tax = integral(f_int_tax,s_min,s_hat);
B_tax = integral(f_int_tax,s_hat,s_star);

int_cost_tax = A_tax + prob*B_tax;

end